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^ ; I. INTRODUCTION 

\ By using direct perturbation theory (ie. strict power series of some coupling constant) to calculate finite tempera- 
ture or nonequilibrium correlation functions, one very frequently encounters infrared (IR) sensitive contributions, ie. 
contributions that become very large (occasionally infinite) when we approach a certain energy scale. These contribu- 
tions, if they show up consequently in every higher loop level, render the direct perturbation theory non-convergent. 
Then we have to collect the problematic diagrams and sum them up, before the actual perturbative procedure begins. 
I , Since these diagrams are usually also ultraviolet (UV) divergent besides their IR sensitivity, the resummation of these 
' diagrams must be accompanied by a proper renormalization. 

In the literature there are several approaches to the renormalization of resummed perturbative series. One of 
• them treats the resummation in the original sense, as sum of diagrams belonging to a certain subset of all Feynman 
dia grams. One possibility is to apply BPHZ renormalization: then, without the reference to the Lagrangian we 
can assign a finite (renormalized) value to every Feynman diagram [1], and we define the renormalized value of the 
sum of diagrams as sum of their BPHZ renormalized value. In principle this method can be applied to all kind of 
resummations, but technically it can be rather involved. The difficulty is to be able to work out the BPHZ forest for 
■>». ' an arbitrary diagram in the relevant subset. This task was performed for the 2PI resummation in various models in 
'■£5 \ [2]- [5] by solving a set of Bethe-Salpeter equations for the four point functions. Diagrammatic analysis can be used 
also to work out counterterms: such a technique was applied in [6] to renormalizc the 2PPI resummation. Diagram 
i| by diagram renormalization can be applied for large N bubble resummation [7], too. 

Another approach is based on the idea that we can try to modify the generating equation of the resummation in 
order to obtain finite results. It is not guaranteed to work, since the equation that we use is, in general, not exact, 
and so we may miss some terms important for UV finiteness. In the Schwinger-Dyson resummation it seems to be 
an important obstacle, cf. [8]. On the other hand, for Hartcc-Fock resummation, such improvement exists, cf. [9]. In 
the O(N) model the resummation is governed by the 1/N order of the diagram. Here we can write up closed forms, 
■ at least in the lowest orders, and so one expects that it can be also renormalized in closed form. In [10] the authors 
present renormalization of the O(N) model up to 0(1/N) order, at the physical point. 

A third approach uses the fact that the contribution of a diagram depends on the renormalization scheme we use. 
It may be possible to find a scheme where the IR sensitive part of all diagrams disappears. In this IR safe scheme 
there is no need for resummation, usual perturbation theory can be used to calculate physical observables. Thereafter 
we can change the renormalization scheme to a standard one (eg. MS) by appropriate change of the renormalized 
parameters of the Lagrangian ( "renormalization group" transformation) . At the end of this procedure we obtain the 
same result in the MS and the IR safe scheme for all quantities up the order (coupling)™ where n is the order of 
calculation. The difference is the resummation effect: it is responsible for the difference that one scheme is IR safe 
while the other is not. 

These thoughts was exploited in [11] to perform a renormalized momentum independent resummation as a change 
of renormalization scheme from an IR safe one to MS. The IR safe scheme was constructed in a way that if in a 
perturbative calculation we encounter IR sensitive contributions, then we choose the finite part of the counterterm 
to cancel this IR sensitive part. This guarantees in a natural way that the calculation is free from IR sensitive 
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contributions. At the end of the calculation we switch back to MS scheme with renormalization scheme transformation, 
thus we obtain the final result in MS scheme. 

This idea is motivated by the "environmentally friendly" renormalization of Ref. [12], where the renormalization 
scale was used to take into account the environmental dependence. 

In [13] the first step was done towards the generalization of the method to momentum dependent resummations. It 
was observed, that if we do perturbation theory with a propagator, which equals to the MS propagator asymptotically, 
then all the overall divergences - and so all the counterterms - are the same as in the MS scheme. Since IR sensitiveness 
occurs - by definition - in the IR, we expect that it can be captured by an appropriate modification in the IR region. 
So if we base our perturbation theory on a propagator that is exact in the IR and MS-likc in the UV, we can use the 
MS counterterms for renormalization. A great advantage of this method is that one does not need to solve Bethe- 
Salpeter equations, which is a central task in other, diagram-based 2PI renormalizations. On the other hand it is not 
easy to implement the momentum independent resummation in this method. 

The purpose of the present paper is to connect the momentum independent and momentum dependent method in a 
common framework. We will examine the possibility and consistency of a perturbation theory with arbitrarily chosen 
free Lagrangian, which not necessarily approaches the usual free part (^(d&) 2 — I ^-$ 2 ) asymptotically. Therefore 
the generated divergences are also different from the usual case. Still it can be consistent, if the usual process of 
the renormalization can be accomplished - the conditions for this will be examined in the paper. But once we 
have a consistent perturbation theory with "any" propagator, then nothing prohibits us to apply a scheme, where 
the renormalized self energy is zero. In this scheme the free propagator is exact - and this is exactly the effect of 
2PI resummation. So we can perform renormalized 2PI resummation with help of a cleverly chosen renormalization 
scheme. 

In this paper we accomplish this strategy for the two-loop 2PI resummation of the <E> 4 theory at finite temperature. 
The paper will be organized as follows. In Section II we study in detail the feasibility of the renormalized perturbation 
theory with an arbitrary propagator, and work out the strategy for 2PI resummation. In Section III we perform 
renormalized 2PI resummation in the quartic scalar model at finite temperature up to two loop level. We close the 
discussion with Conclusions in Section IV. 



II. PERTURBATION THEORY WITH NON-CONVENTIONAL PROPAGATORS 

Let us start from a renormalizable Lagrangian, which now means that it is renormalizable in the non-perturbative 
sense: we can define regularization dependent parameters in the Lagrangian in a way that the exact n-point functions 
are all finite when we remove the regularization. We assume that we split this bare Lagrangian into two pieces: a 
"free" part and an "interaction part" . The free part should be quadratic in the fields in order to be solvable, but there 
are also quadratic pieces in the interaction Lagrangian (like the mass counterterms). Here, however, we do not fix the 
kernel of the quadratic part to the form appearing in the bare Lagrangian: we allow arbitrary choice of the kernel. We 
subdivide the interaction part, as usual, assigning "orders" to them according to a formal loop parameter. In this way 
we defined a perturbation theory. The question is: when is the so-defined perturbation theory finite order by order? 
The counterterms should be defined as the overall divergence of the 1PI diagrams, when "overall divergence" means 
that all the momenta at all lines are going to the infinity at the same rate. When are these counterterms enough to 
make finite the complete contribution of all the diagrams at a given order of perturbation theory? 

For a consistent renormalization, in fact, we should ensure only two conditions [1]. First, the divergence of an overall 
divergent diagram should be local: in Fourier space it should be a polynomial of its external momenta. Second, if 
we subtract all subdivergences of a diagram, the diagram should be overall divergent or finite. The proof that these 
conditions are satisfied in case of using the usual free propagator was done by Weinberg [14]. But the basic property 
of the propagator, which was in fact used by the proof, is that the free propagator can be power expanded "around 
infinity" , ie.: 

G( P + k) = G(p) + k»d^G( P ) + h^d^Gip) + ■■■, (1) 

if p^ — > oo (which means that all of its components go to infinity with a fixed ratio) . This ensures that differentiation 
with respect to external momenta will decrease the degree of divergence of the diagram, and so - if we differentiate 
enough times - the contribution will be finite. The asymptotic expansion property of eq. (1), however, is true for a 
wide class of functions, not just for l/(p 2 —m 2 ). With any of these functions the Weinberg theorem goes through. Then 
all of the consequences will hold true, too, which can be summarized in a phrase that the theory can be renormalized 
by local counterterms with polynomial momentum dependence. 

The Weinberg theorem in itself does not guarantee that we need a finite number of counterterms, only that the 
counterterm method works. We also need the condition that only those operators are divergent that show up in the 
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bare Lagrangian. In the four dimensional <I> 6 model, for example, we generate a logarithmic divergence for the 8-point 
function at one loop level, which was not included in the original Lagrangian, and so it is not renormalizablc. That 
means that we should accept only those propagators, where there are no more divergent diagrams than in the usual 
I/O? 2 — m2 ) case. If we assume that the propagator is asymptotically G{p) ~ p~ Q , then the one-loop correction to 
the 6-point function in the $ 4 model contains the integral ~ f d 4 pp~ 3a . This is UV convergent if a > 4/3: so the 
propagator should be smaller than p -4 / 3 for asymptotic momenta. 

The exact propagator in four dimensions picks up only logarithmic corrections to the 1/fc 2 behavior, according to 
the Weinberg theorem. Thus it may be a good candidate to be the basic propagator of a perturbation theory. 

Once we have fixed the kernel, and so the propagator, too, we can use it for ordinary perturbation theory. In 
particular, we can do finite temperature calculations. Since temperature effects are exponentially suppressed in the 
propagator at high momenta, the overall divergence of a diagram (when all the momenta goes in the same rate to 
infinity) is temperature independent. This ensures that the countcrterms are temperature independent, and so a zero 
temperature renormalization is enough for a finite temperature calculation, too. 

According to the above analysis, a large class of functions is good to be used as "free" propagators. But which one 
is worth to be used? The best choice is the one which is optimized for the problem at hand, ie. where the perturbative 
series converge the fastest. If we keep the bare Lagrangian intact, then any choice describes the same physics, so each 
is a resummed version of another. This is similar when we choose the best renormalization scale. To compute the 
potential, for example, the best choice is to have a scale fi ~ l/r where r is the distance between the point charges. 
Then large logarithmic corrections are suppressed, perturbation series converge well. To keep the bare Lagrangian 
intact we should tune the parameters of the perturbation theory accordingly: in particular we have to work with 
the /i-dependent coupling g(l/r). We will use the same strategy: we optimize the kernel for a certain problem, for 
example for computations at temperature T . To keep the bare Lagrangian intact we should work with T dependent 
parameters in the perturbation theory (g(To)). Then we describe the same physics from different points of view. 

Next we express this strategy in more formal terms in the $ 4 model. First we establish the generic strategy, then 
we examine the relation of perturbation theories with different kernel choices to each other. 



A. Generic strategy 



The bare Lagrangian of the theory reads 

r — ^ (Pl<h \ 2 m barc ,^2 -^barc ^4 (r> \ 
L ~2^ h&Ie > ~ barC ^j~ $ barc- ( 2 ) 

In the usual procedure we first change to the renormalized field: ^arc = Z^ 2 , then we find 

where 5Z = 1 — Z, m 2 , = Zm\ &TC and Ao = Z 2 Abarc- Then we split m 2 , = m 2 R + Sm 2 and Ao = Xr + <5A and find 

£=- <9$ 2 ^$ 2 -<f> 4 -\ (<9$ 2 $ 2 $ 4 . (4 

2 V ' 2 24 2 1 ' 2 24 y ' 

The above form suggests a separation of the non-interacting and interacting theory as 

Co = V) 2 - ^d> 2 , Cr = + -m 2 — $ 2 - -$ 4 , (5) 

u 2 ' 2 24 2 v ; 2 24 w 

where we treat Ci as interactions. It must be emphasized that the values of the counterterms are not fixed beforehand: 
they are determined order-by order, and they do depend on the choice of the chosen perturbative scheme. We can 
prove [1], however, that in the n-th order of perturbation theory the values of the counterterms are scheme independent 
up to 0{X n ). 

The above described procedure is not unique. Starting from (3) we can separate the renormalized and counterterms 

as 

c = hs> K{id)$ - + i $ SK{id) Q g $ 4 (6) 

The kernel of the quadratic part (shortly "kernel" in the followings) used here is arbitrary, although for consistency 
we will require certain properties to be satisfied. To obtain the bare Lagrangian (3), the kernel- and the coupling- 
counterterms should satisfy 

Zp 2 -m 2 = K(p) + SK(p), X = X R + 5X. (7) 



4 



The above form suggests the following free - interaction separation 

Co = l*K*, Cl = -^ + ^SK<P-^. (8) 

The coupling constant and kernel counterterms are considered to be of higher loop order as compared to the coupling 
constant and kernel, respectively. Their values should be determined order by order, just like in the usual case. We 
should compute the diagrams, determine the infinite parts and then make them disappear by a suitable choice of the 
infinite part of the counterterms. The propagator should satisfy (1), then the divergences are polynomial in their 
momenta; it also has to have appropriate asymptotics, then only those operators are divergent that are present also 
in the bare Lagrangian. To make these divergences disappear, we should choose SK^iv = Ap 2 — B and SX — C, where 
A, B and C are divergent constants. Since from (7) we can write S K = SZp 2 — 8m 2 , so we can identify the infinite 
parts of SZ and Sm 2 as A and B, respectively. It must be emphasized, that the actual form of SZ, Sm 2 and <5A depends 
on the choice of the kernel. The logics of their determination is the same as in the usual case when we use the kernel 
p 2 - to 2 . 

B. Theories with different kernels 

The best strategy for perturbation theory would be to use that kernel that is the most appropriate for the description 
of the given observable. The background of working with different kernels is very similar to using different schemes in 
the usual perturbation theory. If the bare Lagrangian is the same, they should describe the same physics: now this 
is ensured by eq. (7). Here the momentum dependence of K and SK should cancel each other, while Z and m 2 are 
determined by the local, divergent part of the kernel counterterm. Similarly, the coupling constant counterterms is 
determined by the divergences of the four point function. 

Although we determine the counterterms order by order, the following statement ensures that in an nth order 
calculation the bare parameters of the Lagrangian, as well as any observables are independent on the choice of the 
kernel, up to 0(X n ). As a consequence, at infinite order, we recover the same bare Lagrangian, and obtain the same 
value for all observables, with any choice of the kernel. Although the proof is essentially the same as in the standard 
case (cf. [I], in particular the Chapter "Oversubtractions" ) it is worth to state here for our somewhat more generic 
case. 

Statement In any regularized perturbation theory with kernel K and kernel counterterm SK which satisfy (7) with 
m o = m 2 R + 5m 2 , Z = l + SZ and Xo = X R + SX where m 2 R and Xr are fixed, and 5m 2 , 5Z and SX are determined 
order by order, the following statements are true in n-th order of the calculation: 

• Sm 2 , SZ and SX are independent on the kernel choice up to order A™ 

• any Greens function is kernel independent up to order X n . 

proof First we observe that if we have a well defined renormalized theory with K and SK, then in a theory with 
K' = K — SK and SK' = SK + SK the bare Lagrangian is the same, and so, nonperturbatively, we obtain the 
same results for all observables. In this case (7) is automatically satisfied for K 1 and SK 1 if it was true for K and 
SK. This proves that in an infinite order perturbation theory the result is independent on the actual choice of 
the kernel. In finite order perturbation theory, however, order-by-order finiteness requires different counterterm 
structure for each kernels. We should prove that the results are stabilized below the order of the calculation 
0(X n ). 

At tree level there is no need for counterterms (SK — 0, SX = 0), and (7) ensures that K = p 2 — m 2 R = K . 

At one loop level we have a one- loop self-energy Si [if], which is only overall divergent. When the regularization 
is fixed then we can determine the "infinite part" of a diagram, so we write 

Si [if] = Sf[if]+Sf [K]. (9) 

In order to ensure finiteness, we should choose the one loop kernel counterterm as 5K\ = Sf lv [if] + 5K\, where 
SK\ is finite (but, possibly, momentum dependent). Then (7) tells us that that the one-loop kernel K\ has to 
satisfy 

ifi + Sf v [ifi] + SK i =K + SZ lP 2 - Sm\. (10) 

We then identify SZ lP 2 - Sm 2 = Sf v [K] and K 1 = K - SK 1 . This latter equation states that K l = K + 0(A) , 
therefore SZ\p 2 — Sm\ — Sf lv [if ] + 0(X 2 ): this means that up to order A it is independent of the choice of 
SK i , ie. of the choice of the kernel. 
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What happens at higher order? We write 8K\ = T,f lv [K] + SKi. This means that we have a new quadratic 
counterterm to the Lagrangian with kernel 5K\. At higher orders this countcrtcrm should be included in any 
diagram T[K). Technically this means that we subsequently replace some propagators by iG — > iG ■ iSK\ ■ iG. 
Since G = K^ 1 , this replacement means 

lG ^^SK ll (11) 
in functional manner. So, up to n inclusions into the diagram, it can be formulated as 

1 d l Y 



E 



l\ dK l 
=o 



{8 Kl f. (12) 



The zeroth term is Y itself, the l/£\ is the remnant of the 1/n! factor of the nth order in the perturbative 
series. We can write it as an infinite sum minus the difference. But the infinite sum is just the Taylor-series of 
Y[K + SK 1 ]= Y[K }. So we find: 

^Timr*®* 1 ?- £ ^|^(^i) £ = r[if + «r 1 ](i-o(A«+ 1 )) = r[x ](i-o(A"+ 1 )). (is) 

i=0 ' £=n+l 

We see, therefore, that the generated finite counterterm restores the value Y[K ] up to order A". This means that 
any renormalized observable, which is represented by the sum of diagrams Y, is the same (kernel-independent) 
up to 0(X n ). Since we did not use any specific property of 5K\ which would exploit its one loop nature, the 
above statement can directly be generalized to higher orders. Moreover, the same line of thought can be followed 
for the coupling constant, when it is momentum independent. This completes the proof of the second part of 
the statement. 

The same proof also applies to the counterterm structure. Then we should use Y — ► £f IV [iT[. At one loop 
level the value of T,f lv [K] defines the kernel counterterms. At higher (nth) loop level the inclusion of 8K\ 
finite counterterm will generate divergent contribution, which has to be made vanish by the nth order kernel 
counterterm. So SK n contains the following contribution 

^-^^™(^r. d4) 

In the sum Y^ii=i therefore, we recover the subseries: 



E 



=o 



{SK^ = ^ v [K } + O(X n+1 ). (15) 



So, although the one loop counterterm depends on the kernel, the higher order counterterm contributions coming 
from SKi inclusions ensure kernel independence up to 0(X n ). Since this line of thought did not depend that 
we start at one loop level, it is valid for any other finite counterterm, and also the same line of thought applies 
to the coupling constant renormalization. This completes the proof of the first part of the statement. 

QED 

In practice, the most important question is the comparability of two perturbation theories, without a reference to 
the bare theory. The above statement claims that this can be ensured if we choose the parameters of scheme "2" 
as m,2(mi, Z\, X\), Z^m-i, Z±, Ai), A2(mi, Z\, X\), where mi, Z\ and Ai are the parameters of scheme "1". These 
relations can be found either by the requirement that the bare Lagrangian is the same (first part of the above 
Statement), or requiring that some observables has the same value in the two schemes. This latter strategy is the 
matching, which is a widely used technique in field theory (cf. large mass decoupling [15] or dimensional reduction 
[16]). It works in our case in a way that we compute three observables from theory "1", and find those m,2, Z%, A2 
parameters that give the same result for these observables from a calculation in theory "2" . The above statement then 
ensures that calculating any other observables, the difference between the results of the two schemes is of 0(X n+1 ) 
where n is the order of the computation. The best is of course to choose that type of observables which can be 
calculated quite reliably from both schemes of perturbation theory. 

The area where we want to apply the perturbation theory with modified kernel is to optimize perturbation theory 
for finite temperature or in non-equilibrium calculations. At each temperatures we can use a kernel, which is optimized 
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to that temperature. To compare the results at different temperatures, however, we have to match the renormalizcd 
parameters. A good set of quantities is the asymptotically large momentum regime observables at T = 0, since we 
expect that these observables should be independent of the details of the kernel at finite momenta, but still sensitive to 
the values of the renormalized parameters. In practical numerical calculations, of course, we should choose quantities 
at large, but not asymptotic momentum, and determine the matching procedure there. 

C. 2PI resummation 

For optimizing the kernel for a given environment, we can use the condition that the free propagator (the inverse 
of the kernel) is the exact one. Then at any higher order diagrams the self-energy insertions, which would modify 
the free - in this case the exact - propagator, should be missing. So we should choose the kernel countcrterms 5K 
at each perturbative order in a way that it exactly cancels the corresponding self-energy T,{K](p), calculated with the 
kernel K. On the other hand the diagrams with self-energy insertions are exactly the 2-particle reducible diagrams. 
These are missing in this perturbation theory: therefore we realized 2PI resummation, differently from the technique 
of Cornwall, Jackiw and Tomboulis [17]. 

In formula, the total self-energy can be written, at nth order, as 

(p) = ^ n) [K](p) 5K^\p), (16) 

where here (and in the followings) E[.K](p) means the self energy with subdivergences subtracted. So 2PI resummation 
requires S^[if](p) = for all n: in this case the exact propagator G~^{p) = K{p) — Y, to t[K](p) = K(p) which is the 
inverse free propagator. So we should choose 

5K^{p) = Y,^[K](p). (17) 

Equation (7) implies 

Zp 2 -m 2 = K(p) + X[K](p) (18) 

at each perturbative order. It seems to be the usual 2PI equation, but it is not quite so. In this formula the infinite 
part of Z and m§ is chosen in a way that it cancels the infinities of T,[K](p): therefore, implicitly, Z and m are 
functionals of K. By our Statement of the previous subsection we can claim that at infinite order this dependence 
vanishes; but at any finite order it is there. 

To see a finite equation we use the fact that E[.K](p) = A^wp 2 — £?div + Sfi n [.K](p), since E, by definition, is just 
overall divergent (subdivergences are subtracted). Then we choose Z = 1 + ^4div + <^C an d m o = m R + ^div, and we 
arrive at (£ = 1 + 5Q 

(p 2 -m 2 R = K(p) + V {in [K](p), (19) 

which is a finite equation. 

We will denote the kernel optimized by 2PI equations for calculation at temperature T as K^ T \ We also will denote 
a diagram or observable V computed with kernel K and at temperature T as T\K, T]. We remark here that we can 
use a kernel, optimized to a certain temperature, in computations at other temperatures as well. The optimization 
temperature is just a notation to label the kernel, while there can be another real temperature in the system. Weinberg 
theorem ensures that r dlv [K , the overall divergent part does not depend on the (real) temperature. 

The strategy to solve the finite temperature 2PI equations will consist of two steps. 

Step 1. We determine the kernel K^> which is optimized for zero temperature, ie. the solution of (19) at T = 0. The 
parameters (, m 2 R and Xr are determined by using physical observables. In the calculations below we just choose 
them to be m R = m, Q = I and A,r = A, so we should solve: 

K {0) (p) =p 2 -m 2 - Sfin^W, 0](p), (20) 

We will use the normalization conditions for the self-energy as: 



(21) 



Then at p 2 = m 2 we will find a pole of the propagator, with unit residuum. 
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Step 2. We fix the temperature at T, and solve (19) to have K^ T \ To be consistent with the T = calculation we 
cannot just choose the renormalized parameters ( m R and \ R ) of the Lagrangian, we will find their value 
by matching to the T = case. 

There are two possible strategies here. The first is to find the solution by taking arbitrary finite part for E, but 
using undetermined values for ( R , m R and X R . Then, using 

K (T) 

as a kernel, we compute the zero temperature 
propagator and 4-point function. Since the kernel was not optimized for T = calculation, we do not expect 
especially good convergence of perturbative series, but we do expect that the asymptotic regime will be well 
described. So we will choose large momenta (eg. p = (Nm, 0, 0, 0) where N = 20 - 50) and match G _1 (p) = K(p) 
and T^(p, 0, — p, 0) at these momenta (choosing together 3 matching conditions). 

The other possible strategy - which will be followed later - defines the finite part of E in a way that it is possible 
to choose Cii = 1 an d m R = m. Then the relevant equation here reads (cf. (19)) 

K (T \p) = p 2 -m 2 - E fin [tf( T ),T](p). (22) 

For matching we still use G^ 1 (p) = K(p) and T^(p, 0, — p, 0) at asymptotic momenta. The zero temperature 
propagator, calculated with the finite temperature kernel reads 

G-^K^Mp) = KV\p) - i:[K {T) A(p) + 5KM- (23) 
Since 2PI equation (ic. E tot = 0) yields 5K x {p) = E[if < T \ T](p), we obtain 

G-^K^Mp) = K^(p) - E[X( T >,0](p) + nK (T \T](p). (24) 

The infinite part of E[K^ T \ 0}(p) and E[if ( T \ T\{p) is the same (overall divergence is temperature independent), 
so we can write 

G- 1 [K^ , 0] (p) = K ^ (p) - E fin [K^ , 0] (p) + E fin [K^ , T] (p) . (25) 
Now we use K {T \p) + E fin [if (T) , T](p) = p 2 - m 2 , then we find 

G- 1 [K^ T \Q](p)=p 2 -m 2 - E Rn [K^ T \0](p). (26) 

This should be equal to 

G-^K^Mp) = Jf(°)(p) =p 2 - m 2 - E fin [^°),0](p). (27) 

So the asymptotic equality of the two propagators requires that we should define the finite part of the self energy 
in a way that 

E fin [^°),0](p) = E fin [^ T ),0](p) (28) 
shall be true for two asymptotic momenta. 
After obtaining the optimal kernel, we can use ordinary perturbation theory to compute other observables. 



III. $ 4 THEORY, TWO LOOP LEVEL, FINITE TEMPERATURE 

Now let us test the above ideas in the case of the $ 4 theory at finite temperature at two loop level. We will use 
real time formalism in R/A basis [18]. Here the original fields (4>i, $2) are replaced by ($ r ,<J> a ) via the definition 

*l=*r + ^*a, $2 = $r-^$a. (29) 

The propagator is a 2 x 2 matrix, where, in this basis G aa = for the exact Green's function, and also G* a (p) = G ar (p) 
is true in the Fourier space. The interaction Lagrangian reads 

-Cj = A$3^ + A$ r $3. (30) 
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A generic kernel is a 2 x 2 matrix, but it must respect symmetries valid for any finite temperature field theory In 
particular, all the propagators should be derived from the spectral function: 

Gn = G ra + G12, G22 = -G ra + G21, iG 12 (p) = n(p )g(p), iG 2 i(p) = (1 + n(p ))g(p), (31) 

and 

00 

G ra (p)= QiUJ,P }. , Q {p) = -2\mG ra {p). (32) 

J 2ir p - lo + is 

—00 

The retarded self-energy G ra can be expressed through its own self energy, it will not mix with other components: 

G; a 1 (p) = G; a \ ( P )-x ar (p). (33) 

The other self energies form a matrix, and they can be obtained as 

Re En = - Re £22 = R- eS ar , ImSu = ImE 2 2 = (1 + 2ns)ImE ar , 

E 12 = -2m B Im£ ar , £ 21 = -2i(l + n B ) ImE ar . (34) 

These relations are to be respected when we choose the generic kernel K. We denote the retarded kernel by K ra , the 
other matrix elements come from (34). The free retarded propagator then simply reads as 

G^oW = Krai?)- (35) 

The free spectral function can be expressed as 



(ReK ra (p)) 2 + (lmK ra (p)) 2 



Causality requires that lmK ra (p) must be an odd function of po, and also lmK ra (p) < must be true for p > 0. 
In the MS scheme it is satisfied by choosing lmK ra (p) = — esgnp (Landau prescription). 

A. One loop level 

At one loop we do not generate momentum dependence for the kernel. This problem has been discussed with 
the momentum- independent version of the present method in [11], but we include it here for completeness and for 
elucidate the strategy. The self-energy and the one-loop four point function read: 

£« [K] (p) = ±T[K] , T< 4 ) (p, q, k, I) = A + y (p + q) + I[K] (p + k)+ I[K] (p + l)} + SX, . (37) 

where T[K] is the tadpole, I[K] is the bubble diagram, computed with kernel K. The tadpole can be expressed as 
Gn(x = 0) = Gi2(a; = 0); in Fourier space 

T\K] = / n{ko)Q(k) = J -0^(1 + 2n(ko))Q(k). (38) 

k >a 

To compute the bubble diagram one should consult Appendix A 1. For the 2PI summation we choose 

8K l = ±T[K]. (39) 

For Step 1 we have to define its finite part so that the self-energy shall be zero on the mass shell at zero temperature. 
It results in the definition 

T[K]=T div [K]+T fin [K}, T div {K}= J j^Q{k), T Rn [K}= J ^2n(k a )g(k). (40) 

fe >0 k >0 
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With this definition the zero temperature equation to solve is (cf. (20)) 

K^(p)=p 2 -m 2 , 

since the self energy is zero at zero temperature. 

For Step 2, with proper finite part for E, we have to solve (22): 



(41) 



(42) 



In order to be write this equation, we have to ensure that the self energies at zero temperature, calculated with 
kernels and are equal. But now the zero temperature self energy, according to (40), is simply zero, so this 

requirement is trivially satisfied. 

(42) is the same result that was obtained also in [11]. 

For the coupling constant renormalization we choose for kernel K the counterterm 



Q\2 

5\ 1 = -—I[K,0](p = 0), 

and define the finite part of the bubble diagram as 

7 fin [K, T] (p, q,kj) = I[K, T] (p, q,k,t)- I[K, 0] (0). 



(43) 



(44) 



This is of course just a specific choice, we could add a finite term to it in another renormalization scheme. The above 
choice ensures that the renormalized one loop four point function at zero momentum at zero temperature is A^ = A. 
We use this equation to write the second matching condition as 



~2 



/fin [K® ,0]( P + q)+ /fin [K (()) , 0] (p + k) + I nn [K™ , 0] (p 



\2 r 



A T + f 



/fin [X (T) , 0] (p + q) + /fi„ [K {T) , 0] (p + k) + / nn [K ^ , 0] (p + I) 



(45) 



at asymptotic momenta. In this simple case only the mass is modified, let us denote the mass of the kernel K as M. 
Then we can explicitly write up these contributions: 



W](0)_ 
/fi„[/C,0](p 2 >4M 2 ) 



1 , M 2 
+ IE + In j 



1 



16tt 2 



In 



P 



where 



X± = 



M 2 



1 ± 



2X + (lnX+ - 1) + 2X_(lnX_ - 1) 



4M 2 



(46) 



(47) 



For asymptotically large momenta (p 2 3> M 2 ) we have 

I Rn [K,T = 0](p) = 



16tt 2 



lnp 2 + ©(const.) 



Therefore, in leading order in momenta, equation (45) can be written (choosing p = q = —k = —£) as 



32tt 2 



lnp 2 = A fl + 



32^ 2 



lnp 2 



(48) 



(49) 



which implies Xt = A. 

We will use this matching condition also at two loop level for simplicity. But we can also argue that - since the 
asymptotic parts of the kernels and are the same (cf. (28)) - we obtain the same /c-dependence for the 

4-point function at asymptotic momenta using either or Then \ R = \ should be the result at all order. 
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B. Two loop analysis 



Let us repeat the above analysis at two loop level. From the generic strategy we assume that the kernel is fixed 
throughout the calculation. The diagrams for the self energy read: 

^[K](p) = ^-S[K](p) + ^T[K] I[K}(0) ^I[K](0) + S -^T[K], (50) 

where T[K] stands for tadpole, I[K] for bubble, S[K] for setting sun diagrams. The computation of the setting sun 
diagram is discussed in Appendix A 2. 5K\ and 5Xi is fixed at one-loop level (cf. (39) and (43)). 

To see the overall divergence we split the setting sun diagram as S[K](p) = Ss n [K] (p) + S su bdiv[K] + S oveia \\[K](p) , 
where 

S suhdiv [K] = 3T[K]I[K, 0](0) = -25MTIK], S oveiM [K](p) = SZ sP 2 - dm 2 , (51) 

and 5Zs and 5ms depend on the kernel and the choice of the renormalization conditions. We find 

\2 n 
^ 2) [K](P) = ^rS^[K]{p) + 5Z sP 2 - 5m 2 + -±T[K], (52) 

6 

We want to make the total self-energy correction disappear in the spirit of the 2PI resummation (cf. (17)), so we 
should choose 5K2(p) = 'E^[K](p). To tell the bare parameters of the Lagrangian we should write (cf. eq. (7)) 

Zp 2 - ml = K(p) + 5K 1 +5K 2 . (53) 

We write mf t = m\ + 5m\ + 5m\ and Z = 1 + 5Z\ + 5Z 2 . Then we identify 

\ A \ 

6ml = = -7T T div, 6ml = 6ml - -±T[K], SZ 1= 0, SZ 2 = -6Z S . (54) 

2 6 

To renormalize the theory we should perform the two steps described earlier: 
Step 1 we have to solve the zero temperature 2PI equation (20): 

#<°>(p) =p 2 -m 2 ^ fin [^(°),0](p), (55) 



since Ta n [K^°\0] = 0. We should define the finite part as 

S Rn [K^,0](p 2 = m 2 ) = 0, 



dp 2 



0. (56) 



Then the above equation yields a propagator which has a pole at p 2 = m 2 with unit residuum. 

In numerical calculations one has to compute the complete contribution, subtract 3T[K^°\0]1[K^°\0\(0) subdi- 
vergence, and subtract an ap 2 — b function. The value of a and 6 can be determined from eq. (56). 

Step 2 the finite temperature equation to solve is (22): 

K^{p)=p 2 -m 2 - ^[K^\T] ^ Siin [K( T \T](p). (57) 

For the correct renormalization we should require (28): 

S Rn [K (T) , 0](p) = S Rn [K^\ 0](p) for asymptotic momenta. (58) 

Numerically we should compute the complete setting sun contribution with kernel at zero temperature, 

subtract ?>T[K {T \ 0]J[K^ T \ 0](0) subdiver gence, and finally Ap 2 — B. We determine A and B from (58), by requiring 
its fulfillment at two asymptotic momenta pi and p 2 ■ Then at finite temperature we should compute the rcnormalizcd 
value of the setting sun diagram, that we first compute the complete diagram with kernel at finite temperature, 
subtract 3T[K( t \T]1[K<- t \0}(0) subdiver gence, and also Ap 2 — B with the values determined before. This will cancel 
all divergences, since, as it was mentioned earlier, the overall divergences do not depend on the temperature. 
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C. Results 



The above strategy can be accomplished by numerical calculations. We used the spatial rotational invariance of 
the propagators, which allows to represent them at finite temperature as two-dimensional functions, with variables 
Po and |p|. Therefore for the integrations we need a 2D lattice, which is mapped with a continuous function to the 
infinite 2D space. All the results presented here are calculated at A = 10. Finite temperature, if not stated otherwise, 
means here T = m. 

The 2PI equations are solved by successive approximation, starting from a spectral function near to the free one. 
From the spectral function one can determine the propagators, with the propagators we calculate the self-energies, 
and with help of self energies we update the spectral function. Then repeat this algorithm till it converges. 

First we discuss the renormalization. In Fig. 1 we can see the figures representing the renormalization process. On 
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FIG. 1: Quality of renormalization: a.) zero temperature self-energy at p = 0; b.) zero temperature setting sun diagram with 
different kernels - the two lines are almost within line width. 

Fig. 1/a there appears the self energy calculated with . It can clearly be seen that at po = m the self energy is 
£ = 0, so there is a pole there. It can also be seen that the derivative with respect to po is zero - this ensures the unit 
residuum. Fig. 1/b concerns the finite temperature renormalization. Here we require that the values of setting sun 
diagram calculated with kernel or at zero temperature, are, asymptotically, the same. As this plot shows, 

this is fulfilled within line width, although for smaller momenta there can be seen some discrepancy. These two plots 
show that the theory is renormalized. 

The next two plots (Fig. 2) show the spectral function at two different spatial momentum. At zero temperature one 




can clearly see the quasiparticle peak, and a well separated threshold starting at 3m in case of p = 0. The position 
of the peak as well as the threshold value shifts with the spatial momentum. At finite temperature the quasiparticle 
peak is still very sharp and narrow (at these values of the coupling), but there is no sign for any threshold behavior, 
only a broad continuum can be observed. 
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In the third pair of figures (Fig. 3) the quasiparticle properties can be seen. The first plot shows the spatial 
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FIG. 3: The quasiparticle properties: spatial momentum dependence of the mass (a.) and the width (b.). 

momentum dependence of the mass. The momentum dependence is rather weak: in the range < |p| < 15m it is 
less than 10%. This fact, together with the dominant quasiparticle peak may explain the success of the "screened 
perturbation theory" [19]. In the second plot the quasiparticle width can be seen: it drops[21] about a factor of 4 in 
the range < |p| < 10m, as was predicted by perturbative calculations [20]. 

The last two plots demonstrate the temperature dependence of the quasiparticle properties. In Fig. 4/a. the 
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FIG. 4: The temperature dependence of the quasiparticle properties: mass (a.) and the width (b.). 



temperature dependence of the on-shell mass can be seen. In the range T < 0.5m there is hardly any temperature 
dependence; at T > 1.2m the curve fits well to a second order polynomial, so there m 2 = m 2 , + cT 2 is a good 
assumption. Fig. 4/b shows the on-shell damping rate as a function of the temperature. The same tendency can be 
observed as in the on-shell mass case: for small temperatures (T < 0.5m) there is hardly any damping can be seen, 
while for high temperatures (T > m) the damping rate can be fitted to a straight line T = r + cT. In both cases the 
T « m regime is a crossover between the two characteristic behavior. 



IV. CONCLUSION 



The most important message of this paper is that the "renormalization scheme method" is appropriate to perform 
fully renormalized momentum dependent 2PI resummation at finite temperature. The advantage of this method is 
that one does not need to solve the Bethc-Salpctcr equations. Instead one fixes the renormalization in the finite 
temperature calculation by matching to the zero temperature results. 

Renormalization scheme method is based on the observation that in certain momentum dependent schemes there 
is no need of doing 2PI resummation, since, by construction, the self energy is zero at each order. We studied in the 
paper, when such an "optimized" scheme yields a consistent perturbation theory. We examined the relation of two 
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such schemes, and proved that by appropriate choice of the renormalized parameters, two schemes can give the same 
result up to 0(X n ), where n is the order of the perturbation theory - this is the same result as in the "ordinary" 
momentum independent schemes. But it also suggests that if we fix a scheme at zero temperature, then at finite 
temperature the optimized scheme can be related to it by matching. 

In particular we suggested in the paper a two-step method for the renormalization. In the first step one has to solve 
the zero temperature 2PI scheme equations, where we fix its parameters by observable quantities (like quasiparticlc 
mass, residuum and scattering cross section). In the second step we solve the finite temperature 2PI scheme equations, 
defining the finite part of the self energy in a specific way: we require that the zero temperature self energy, calculated 
with the zero temperature and the finite temperature optimized kernel, are the same for asymptotic momenta. This 
definition is the manifestation of the matching procedure. 

The so-defined rules are numerically easily implementable, as it was demonstrated in the paper on the example 
of the <I> 4 theory with two-loop 2PI resummation. The results of this study are summarized in the paper. The 
most important corollary is that the $ 4 theory can be well characterized by a thermal mass: the spectral function is 
dominated by a quasiparticle peak where its spatial momentum dependence can be described by an almost constant 
mass term. 

For future studies, we plan to generalize the renormalization scheme method to non-equilibrium case, and also - 
which may be even more important and interesting - to systems with specific symmetries (global or gauge) . 
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APPENDIX A: DIAGRAMS 



1. 



Bubble 



We compute the retarded diagram that reads 




(Al) 



Using the dispersion relation 




(A2) 



the retarded diagram can be computed from its discontinuity 




(A3) 



At zero temperature 



n(ui) = 



1 



—> -e(-w). 



(A4) 



Thus 



1 + n(p ) + n(k a - p Q ) = 0(p o ) - ®(po - fco) = ©(0 < Po < k ). 



(A5) 



So 




(A6) 
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2. Setting sun diagram 



The setting sun diagram can be reproduced from its discontinuity, just like the bubble. The discontinuity reads 
meiS[K\(k)= J-^-^[iG 21 (p)iG2i(q)iG 21 (e)-iG 12 (p)iG 12 (q)iG 12 (e)) 

[(l + n p )(l+n q )(l + ne)-n p n q ne]Qi(p)Q2{q)Q3{t), (A7) 



(2tt) 4 (2tt) 4 

where p + q + £ = k. At zero temperatures it falls back to 



DisczS[K,T = 0](fc) = | ^jL g 1 (p)g 2 (q)g 3 (e) [O(p )O(q )Q(e ) - Q(-p )Q(- qo )Q(-e )] . (A8) 
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